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The effect of dipolar and exchange interactions within pairs of paramagnetic electronic states on 
Pauli-blockade-controlled spin-dependent transport and recombination rates during magnetic reso- 
nant spin excitation is studied numerically using the superoperator Liouville-space formalism. The 
simulations reveal that spin-Rabi nutation induced by magnetic resonance can control transition 
rates which can be observed experimentally by pulsed electrically (pEDMR) and pulsed optically 
(pODMR) detected magnetic resonance spectroscopies. When the dipolar coupling exceeds the 
difference of the pair partners' Zeeman energies, several nutation frequency components can be ob- 
served, the most pronounced at ^/2jBi (7 is the gyromagnetic ratio, Bi is the excitation field). 
Exchange coupling does not significantly affect this nutation component; however, it does strongly 
influence a low-frequency component < jBi. Thus, pEDMR/pODMR allow the simultaneous iden- 
tification of exchange and dipolar interaction strengths. 



I. INTRODUCTION 

In solids with weak spin-orbit coupling like silicon or 
carbon-based materials, spin-selection rules induced by 
spin conservation can drastically influence optical and 
electrical materials properties^ ^. Because of this, a ma- 
nipulation of spin states, e.g., by means of magnetic res- 
onance, can change conductivity, luminescence, or ab- 
sorption. These effects can be used for the investigation 
of the microscopic physical nature of the paramagnetic 
species involved in these processes, as is done with ex- 
perimental techniques such as electrically (EDMR) and 
optically (ODMR) detected magnetic resonance spectro- 
scopies. An abundance of spin- dependent processes has 
been reported in the literature^^^^^^^. Most of these re- 
ports involve the Pauli-blockade effect, where a transi- 
tion of two paramagnetic states with s = 1/2 into a 
single doubly occupied electron state with singlet con- 
figuration is controlled by the singlet content of the pair 
before the transition occurs. Such mechanisms are usu- 
ally described by an "intermediate-pair" process, where 
an exclusive pair of two spins is formed. This pair then 
either dissociates with spin-independent probability, or 
undergoes a transition into the singlet state that happens 
with probability ex KS'I^)^, where |^) is the spin state 
of the pair before the transition. This intermediate-pair 
model, developed by Kaplan, Solomon, and Mott (KSM) 
in 1978^, is distinct from other s = 1/2 pair models 
which do not require the exclusive intermediate pair (see 
for instance t he Le pine^ model). However, many experi- 
mental EDMRpE^ and ODMRP studies have shown the 
validity of this picture for the description of several spin- 
dependent recombination and transport effects which in- 
volve transitions between localized electronic states. The 
KSM model is thus the basis for the calculation of spin- 
dependent transition rates presented here. 

With the availability of high-power microwave sources, 
and the resulting development of pulsed electron param- 
agnetic resonance (EPR) techniques in the past 25 years. 



ODMR and EDMR have increasingly been conducted as 
transient, pulsed (p) experiments, on time scales where 
coherent spin- motion effect^^^^^ take place. Since co- 
herent propagation of a quantum mechanical system is 
directly controlled by its Hamiltonian, this development 
in experimental techniques has dramatically enhanced 
access to the fundamental physical nature of the micro- 
scopic systems responsible for the EDMR and ODMR sig- 
nals. Coherent spin effects such as spin echoes, spin-Rabi 
nutation or dynamic decoupling schemes have produced 
a variety of experimental insights. In order to fully utilize 
the spin effects observed with these techniques, a rigorous 
theoretical understanding of the signals is necessary. As 
pEPR spectroscopy evolved over the past decades, many 
studies have contributed to the development of this un- 
derstanding^^ and it is now straightforward to derive 
information from pEPR data about Lande factors, spin- 
spin coupling phenomena, such as exchange coupling, 
dipolar coupling (which reveals distance between inter- 
acting spins), hyperfine couplings, and relaxation times, 
among other variables. Unfortunately, this rather com- 
prehensive theory of pulsed EPR spectroscopy is only 
partially applicable to pulsed EDMR and ODMR exper- 
iments. EDMR and ODMR are performed by measure- 
ment of spin-dependent rates whose observables depend 
on the permutation symmetry of the involved spin pairs, 
i.e., the singlet and triplet content. Thus, the observable 
of EDMR and ODMR experiments is fundamentally dif- 
ferent than the observable of EPR experiments, the latter 
being the magnetic polarization of the observed spin en- 
semble. As a consequence, a spin ensemble that is simul- 
taneously observed with EPR and EDMR/ODMR may 
exhibit entirely different signal behavior due to the dif- 
ferent observables onto which the observed spin ensemble 
is projected. 

Several recent studies aimed at developing and un- 
derstanding pEDMR and pODMR signals have focused 
on electrically or optically detected transient nuta- 
tion measurements, where a spectroscopy of observed 
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spin-Rabi oscillation is conductecP^^^tlll^ These stud- 
ies have considered various spin-coupling regimes for 
the spinpair, including the absence of any spin-spin 
coupling ^ l ^^ l ^^ l, the presence of exchange interaction^^, 
and a disorde r- indu ced distribution of spin-orbit interac- 
tion strength j-^^ l ^^ i Recently, the first analytical study 
of coherently controlled spin-dependent transition rates 
within pairs of strongly exchange- and dipolar-coupled 
pairs was conductecPH. However, a general numerical 
or analytical study for electrically or optically detected 
transient nutation of pairs with arbitrary spin-dipolar 
and spin-exchange interactions is lacking. Such a study 
is the focus of the present work. 

II. INTERMEDIATE- SPIN-PAIR MODEL WITH 
DIPOLAR AND EXCHANGE COUPLING 

Following previous discussions of sp in- depende nt tran- 
sition controlled by intermediate pairJ^HEHIini]^ 
scribe the dipolar- and exchange-coupled intermediate- 
spin- 1/2 pairs with the Hamiltonian 

^spin = h[B • (jaSa^lbSb) - JSa • ^5 

Here, the first term represents the Zeeman interaction for 
both spin-pair partners, the second term is an isotropic 
exchange interaction, the third is a secular (high-field 
approximation) magnetic-dipole coupling, 7^ and 75 are 
the effective gyromagnetic ratios of the spin-pair partners 
a and 6, respectively. The magnetic field 

B=zBo^xBie-''^^ (2) 

consists of a static component Bq along the z-axis, and 
an oscillating component that is chosen to be along the 
a^-axis. When the exchange interaction strength J and 
the dipolar coupling strength D are scaled by they 
can be directly compared to the Larmor separation Au 
of the electron and hole. The negative signs in front of 
the D and J terms are chosen to represent an attractive 
electron-hole paii^^ We note that changing the sign of 
sign of J and/or D will not change the results presented 
below (such a sign change could occur for like-charge spin 
pairs, e.g., bipolarons^). 

The spin-pair Hamiltonian in absence of radiative ex- 
citation {Bi = 0) is rotated into an energy eigenbasis by 
a Jacobi rotation, Hq^ = [/^i^spin^, with the resulting 
eigenbasis given by 

/in)\ / in) \ 

rrt in) _ COs((/>)|n)-sin((/>)|;t) 

^ |;t) - cos((/>)|n)+sin((/>)|;t) ' 

\\u)/ \ \u) J 

where cot(2^) = In the case of either strong dipo- 

lar or strong interaction exchange, the energy eigenbasis 
becomes a set of singlet and triplet states. With strong 



dipolar coupling, (j) ^ —j {D ^ 00, J^O), the en- 
ergy eigenbasis becomes {|T+) , |To) , 16*) , |T_)}; strong 
exchange coupling, (j) ^ j (J^oo, Z^^O), produces 
an energy eigenbasis {|T+) , 15) , |To) , |T_)}. In either 
one of these strong coupling cases, the only ESR- allowed 
transitions are those within the triplet manifold, lead- 
ing to a strong triplet ESR signal. However, because the 
triplet-singlet transition probability is zero, there is no 
observable pODMR/pEDMR signal. Any intermediate 
case (e.g., J ^ D ^ Au) will have an energy eigenbasis 
of {|T+) , |2) , |3) , |T_)}, where |2) and |3) wih each have 
a mixture of singlet and triplet content defined by the rel- 
ative magnitudes of the dipolar and exchange strengths. 
Therefore, the transitions between states are uniquely 
governed by the collection of system parameters J, 
and Auj. 

For pODMR/pEDMR experiments on intermediate- 
spin-pair processes, the observable depends on the per- 
mutation symmetry of the individual pairs, contrary to 
most conventional spectroscopy experiments which probe 
polarization states. An extended discussion of such 
intermediate-pair related pEDMR/pODMR observables 
is given by Gliesche et al.^, who established the con- 
nection of the spin density operator p of an ensemble of 
spin- 1 pairs to a spin-dependent rate transient 

Q(r)= ri?(t)^t = ^(p,,(r)-pfj(l-e-^^^o), (4) 
•^0 i=i 

which follows coherent spin excitation. In Eq. [4j the 
density matrix is in the 4x4 energy eigenbasis repre- 
sentation and the time-dependent function R{t) is the 
spin-dependent rate after the pulse excitation, which 
is assumed to end at t = 0. Since R{t) is a current 
for pEDMR experiments, the integral Q{r) becomes a 
number of charge carriers which undergo spin-dependent 
transitions due to the resonant spin excitation. The de- 
pendence of Q on the pulse length r will reveal infor- 
mation about how the density operator p evolves from 
the steady state to a coherent state due to the resonant 
excitation. Thus, Q{r) is an easily accessible observable 
for the coherently manipulated spin ensemble, represent- 
ing either the number of charge carriers (for pEDMR) or 
photons (for pODMR). 

The transient evolution of Q{r) during the pulse can be 
Fourier transformed (FFT{(5(r)}) in order to make the 
frequency components of the coherent spin motion ex- 
plicit. A comparison of experimentally obtained Rabi fre- 
quency spectra with calculations we present below gives 
insight into the nature of the spin-pair Hamiltonian. As 
the spin-pair Hamiltonian crucially depends on the mi- 
croscopic nature of the spin pairs, pEDMR/pODMR ex- 
periments are superb probes to gain unambiguous exper- 
imental access to spin-dependent transport and recombi- 
nation processes. 

Again following previous descriptionJ^E^t^H^, we de- 
scribe the evolution of the density operator p by a 
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stochastic Liouville equation 



dtp=-[p,H,A+S[p], 



(5) 



where the stochastic term S[p\ =5'cr[p]+'S'an[p] is the sum 
of creation and annihilation terms of the spin pairs. 

As shown elsewhere^, the recombination probabihties 
for the different energy eigenbasis states are given by 
ri = rs \{i\S)\ + tt \{i\T)\ , where rs and rx are the 
singlet and triplet recombination probabilities, respec- 
tively. Using Eq. |3) the various recombination rates can 
be expressed in terms of the coupling parameters by 



1 / J-D \ 

2 I V(J-D)2+(?iAw)2 



(6) 



(7) 



The eigenstates |1) and |4) always remain pure triplet 
states (|T+) and |T— ), respectively); their recombina- 
tion rates are thus not affected by any coupling within 
the spin pair. Under strong coupling such as I) ^ Acj, 
r2 = Vt and rs = rs {J ^ Acj, r2 = Vg and = tt). The 
dissociation rate coefficient d is assumed to be spin inde- 
pendent. In the energy eigenbasis, the stochastic anni- 
hilation term San[p] has matrix elements in a convenient 
form, {San[p]}ij = (^i + + 2d)^. We also assume 
pair generation only creates pairs in an energy eigenstate, 
{^cv[p]}ij = where k is the net generation rate of 

all four states. This creation term is the only inhomoge- 
neous contribution to Eq. [5] In this paper we neglect the 
Redfield relaxation matrix, an assumption that is valid 
in the short-time regime (r < ^ ~ T2 < ^ <C Ti). For 
the purpose of obtaining sufficient resolution, some pulse 
lengths violated this assumption (see Fig. [TJd) and (g)). 



methods make the convolution of distributions of many 
parameters (J, bandwidth of pulse, etc.) impractical 
without a supercomputer. We use two techniques that 
lead to a significant decrease in computation time. In 



Sect. |III A| we detail the first step of the computation, a 
transformation into the rotating frame. 

Once in the rotating frame, several limiting cases 
of the Rabi nutation frequencies are demonstrated in 
Sect. IIIB The limiting cases of overall weak coupling, 
strong exchange coupling, strong dipolar coupling, and 
a large difference in dipolar and exchange coupling are 
described. Sect. [Ill A| also includes an analytical descrip- 
tion of the y/2^Bi Rabi frequency component that oc- 
curs in the presence of strong dipolar coupling. These 
limiting cases provide significant insight into qualitative 
features observed in the numerically calculated general 
gases, such as resonance location, Rabi frequency, and 
signal amplitude. 

In addition to the use of the rotating frame, the calcu- 
lation of the time-dependent change of the density ma- 
trix was aided by the use of Liouville-space formalism, 
and is discussed in Sect. Ill C A direct consequence of 



this formalism is that the inhomogeneous stochastic Li- 
ouville equation is cast into a readily tractable and solv- 
able form. Compared to previous worli^^^^, the speed 
of the simulation allows us to perform a larger and more 
detailed study of Q(r)'s dependence on dipolar D and 
exchange J interactions with respect to the Larmor fre- 
quency separation Acj and the excitation-field strength 

In Sect. |IV[ representative results of these simula- 
tions are given and discussed. Using the methods from 
Sect. Ill, we simulate Q(r) for a range of values D and 
Acj with a fixed excitation field Bi. Then, Q{t) is sim- 
ulated as a function of D and J with a fixed Acj and 
El. Finally, Q(t) is simulated with small and large ex- 
change coupling strengths, along with a complete Pake 
distribution of dipolar interaction strengths. 



III. ANALYTICAL AND NUMERICAL 
METHODS 

In the following section we outline our study of the ob- 
servable Q{t) that results from the coherent excitation 
of the spin pair. Eq. |5] is a set of sixteen coupled inho- 
mogeneous ordinary differential equations (ODEs) that 
were previously solved using a Runge-Kutta or compa- 
rable ODE solvei*^^^^. These computationally intensive 



A. Rotating- Frame Stochastic Liouville Equation 

The rotating frame corresponds to a transformation 
of the Hamiltonian from the energy eigenbasis: Hr = 
R^Hq^R. The rotating-frame density matrix is then 

given by pR = R^ pR. Here R = Ri ^Ri is the 4 x 4 
spin-^ pair rotation operator, and Ri is the rotation op- 
erator for a spin-^ state around the z-axis by an angle 
ujt. The resulting rotating-frame Hamiltonian is 
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hsin(^)) 



V 



2wo — f— -D 7i?i(cos((/)) — sin((/))) 'yBi{cos{(j))+^iiiy<y)) u 

7Bi(cos(<^)-sin(,^)) j + D+y/Auj^ + {J-Df 7Bi(cos(<^)-sin((/))) 
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(8) 
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and has no explicit time dependence. Note that the en- 
ergy levels for the energy eigenbasis Ei, ^3, and E4 
reside on the diagonal. We label the average of the spin- 
pair Larmor frequencies ujq = {uja -\- oJb) and Larmor 
frequency separation Auj = uja — ^b- We assume that the 
Rabi frequency of each spin is the same (7a^i ~755i), 
allowing us to explicate the results in terms of a single 
on-resonance Rabi frequency 7^1. Neglecting the very 
small difference in the individual-spin Rabi frequencies 
symmetrizes the simulations about u — uq = 0, rather 
than demonstrating an inconsequential asymmetry. Af- 
ter an additional time-independent term F = R^dtR is 
absorbed into an effective Hamiltonian H = Hr—F^ the 
rotating-frame stochastic Liouville equation becomes 

dtPR = \[pR.H]^S[pR]. (9) 

As expected from this transformation, the only term left 
with time dependence in Eq. |9]is the rotating-frame den- 
sity matrix pR. 

B. Limiting Cases of the Rabi Frequencies 



transition that has a degeneracy of four and a Rabi fre- 
quency 7^1 . However, if there is a sufficient excitation- 
field strength Bi, both uncoupled spins will nutate co- 
herently, creating a spin-beating effect with a 27^1 Rabi 
frequency component^ 

For the second limiting case consider an effectively 
weak coupling, where the difference in coupling strengths 
becomes much less than the Larmor separation, \J—D \ <C 
Acj. In this limit there are four non-degenerate resonant 
transitions. As in the weak regime, a pair in the effec- 
tively weak regime has a Rabi frequency 7^1 equal to 
that of a single uncoupled spin. 

Both weak and effectively weak coupling leave the en- 
ergy eigenbasis completely unaffected by the rotation 
performed in Eq. [3] In the latter case, this happens 
even though the couplings J and D could individually 
be quite large compared to Acj. However, the resonance 
frequencies for each transition will be shifted due to the 
increased coupling strengths. This non-degenerate en- 
ergy spectrum distinguishes the effectively weak coupling 
from weak coupling; see Fig. [3j 



Useful equations that elucidate limiting cases can be 
derived from finding the single-transition Rabi frequen- 
cies of the rotating-frame Hamiltonian given in Eq. [S] 
By considering an induced transition between only two 
of the available four states and solving a 2 x 2 eigenvalue 
problem, it can be shown that the single-transition Rabi 
frequencies are 



(l^sin' 



(10) 



The negative sign in the first term under the radical on 
the right hand side gives the Rabi frequencies for the 
\T±) ^ 1 2) transitions between the pure triplet states 
and |2) state {{ij) = {(T+, 2); (T_, 2)}). The plus sign 
in Eq. [To] gives the Rabi frequencies for the \T±) ^ |3) 
transitions between the pure triplet states and |3) energy 
eigenstate {{i^j) = {(T+, 3); (T_, 3)}). In general, there 
are four resonant frequencies, ujij = {Ei — Ej)/h. If an 
on-resonant excitation frequency uj is applied such that 
uj = cJij, the second term in Eq. IlO] vanishes. 



1. Weak and Effectively Weak Coupling 

For the first limiting case, let the coupling terms J 
and D approach zero. In this weak-coupling regime, 
I J|+|D| <C Acj, the first term in Eq. 10 tends towards the 
limit (l^sin2^) 



1 and there is only an on-resonance 
Rabi oscillation frequency of a single uncoupled spin, 
7^1. There are two resonant transitions with a two- 
fold degeneracy corresponding to the Larmor frequency 
of each spin in the pair. If the Larmor separation Ao; 
is zero (indicating that the gyromagnetic ratios of the 
electron and the hole are the same), there is only one 



2. Strong Dipolar Coupling 

Now consider the limiting case of strong dipolar cou- 
pling, \D\^ with no exchange coupling, J = 0. As 
D gets large, sin 20^—1, and the four resonant single- 
transition frequencies, offset from cjq, are approximately 



The first term under the radical on the right hand side 
of Eq. [To] is (1 - sin 20) ^ 2 for the |T±) ^ |2) transi- 
tions, and (1 + sin 20) ^0 for the \T±) ^ |3) transitions. 
Therefore, strong dipolar coupling within the pair yields 
an on-resonance Rabi frequency of \/2jBi for each tran- 
sition between the pure triplet states and the |2) state. 
The T±^\2) transition probabilities are large but have 
an overall reduction of the pEDMR/pODMR signal, ow- 
ing to the strong triplet character of the |2) state. 

We predict a Rabi frequency for any spin-^ 

pair with sufficient Larmor separation and strong enough 
dipolar coupling. When strongly coupled, an applied ex- 
citation necessarily affects both spins in a pair, even if 
only a monochromatic excitation is applied. The strong 
dipolar coupling (like a strong exchange coupling^^) al- 
lows access to only one quantum state, and prohibits 
isolating an individual spin within the spin pair. This 
behavior has been well known from traditional magnetic 
resonance spectroscop}'^^ and, without explicit theo- 
retical p roof, it has already been applied to experimental 
pODMlP™ and pEDMR^^ data. 
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3. Strong Exchange Coupling 

We now consider the strong exchange couphng regime, 
where \J\ ^ IAcjI, with no dipolar coupling, D = 0. As J 
gets large, sin 20 1, and the resonant single-transition 
frequencies, offset from cjq, are approximately 



4J 



^±,3 



" 4J ' 



(12) 



The first term under the radical on the right hand side of 
Eq.[lO]is (l+sin2(/))^2 for the |T±)^|3) transitions, and 
(1-sin 20)^0 for the |T±)^|2) transitions. The single- 
transition analysis predicts a Rabi frequency for 
the |T-t) ^ |3) transitions. However, this naive analysis 
does not take into account that the splitting in |T±) ^ |3) 
transition frequencies is so small that the transitions will 
be driven simultaneously by Bi. Therefore, Eq. [To] is 
no longer valid, and a multiple-transition analysis must 
be used. The \T±) ^ |2) transitions are far away from 
ujQ and have Rabi frequencies approaching zero. Because 
of this, the strong exchange-coupling regime can be ana- 
lyzed using only the |T-t) ^ |3) transitions. If a excitation 
frequency oi uj = uoq is applied, two of the three rotating- 
frame energy eigenvalues in the multiple-transition anal- 
ysis are degenerate. This simplifies the eigenvalue prob- 
lem significantly, and the (three-state) Rabi frequency is 
found to be 



•2(l+sin20)(75i)2 ^ 27^1 



(13) 



for the |T+) ^ |3) ^ |T_) transition. We note that the 
single-transition Rabi frequencies for |T-t) ^ |3) do not 
merely add to a total 2V^7^i Rabi frequency. 

If the power of the excitation field is lowered, only 
a 7^1 Rabi frequency is seen in the case of uncoupled 
pairs, whereas a strongly exchange-coupled pair always 
has a 27^1 component, provided the signal is strong 
enough. This fact has served to distinguish uncoupled 
and s trong ly exchange-coupled states in experimental 
studieP^. 



4. Large Difference in Exchange and Dipolar Strengths 

The final limiting case we consider is to take the dif- 
ference in coupling strengths to be large with respect to 
the separation of the Larmor frequencies, and the ex- 
change strength to be greater than the dipolar strength, 
J — D^ Auj. In this limit we have (1 — sin2(/)) for 
the \T±) ^ 1 2) transitions, and (1 + sin 20) 2 for the 
\T±) ^ 1 3) transitions. The resonant single-transition 
frequencies, offset from cjq, are now 



^±,3 



± 



D 

2 ^ 4( J-D) 

3D 

2 4( 



i{J-D) ) ' 



(14) 



The presence of dipolar coupling splits the transition fre- 
quencies enough that the single-transition analysis for 
Eq. [8] becomes valid again. Therefore, in the limit 
of a large difference in dipolar and exchange coupling 
strengths, a Rabi frequency of y/2^Bi will occur when 
on resonance with the \T±) ^ |3) transitions, and the 
|T±) ^ |2) transitions have a vanishingly small transition 
probability. 

The other limit in the strong-coupling regime that we 
do not describe in detail is the difference in coupling 
strengths large with respect to the separation in Larmor 
frequencies, and the dipolar strength is greater than the 
exchange strength, D — J^Acu. An analysis similar to 
that given above shows that the V2jBi Rabi frequency 
components exist, but have resonances far away from the 
central average ujq of the spin-pair Larmor frequencies. 

We will refer to these limiting cases as we discuss the 
features appearing in the results of the following simula- 
tions. 



C. Liouville- Space Formalism 

We now reformulate the rotating-frame description us- 
ing Liouville operator space, also known as superopera- 
tor formalisnP^i^, to increase the computational power 
of the simulation. This technique was also used recently 
in a model for magnetic-field effects in disordered semi- 
conductorJ^il. The essence of this reformulation is the 
representation of the state population as a 16x1 column 
vector p instead of the typical 4x4 density matrix p. Op- 
erations involving A are associated with corresponding 
superoperators A. Note that this formalism produces no 
new physics, but simply recasts the problem such that a 
convenient, tractable solution to Eq. |9]is obtained. 

Using superoperator formalism, the rotating-frame in- 
homogeneous stochastic Liouville equation (Eq. |9| can 
be rewritten in the compact form 

dtpR=^HpR^S^npR + K = GpR^K. (15) 

Here HpR is the abbreviated superoperator form of the 
commutator [pn^H]. is a 16 x 16 superoperator that 
can be written as 



H-- 



( H-IHn 



IH12 
H — IH22 

IH42 



IH2A 

IHs4 

H-IHm J 

(16) 

where / is the 4x4 identity matrix and Hij are the ma- 



IH23 
H — IH33 



\ 



trix elements of the 4x4 Hamiltonian H. In Eq. 15 
5'an is a time-independent diagonal 16 x 16 matrix of the 
appropriate stochastic annihilation terms corresponding 
to Sg,n[p]' The creation term i^T is a time-independent 
16 X 16 matrix consisting of the appropriate stochastic 
creation/generation terms corresponding to Scr[p] and is 



6 



the sole inhomogeneous part of Eq. 15 The superop- 



erator G is merely the addition of and 6'an; it is a 
symmetric and relatively sparse matrix with 160 zeroes. 

A steady-state density matrix is used to define p(0), 
the density matrix at time t = 0, and is obtained by 
neglecting the coherent excitation {Bi = 0) and finding 
a steady-state superoperator Gs from Eq. [TsJ Using the 
variation-of-parameters method, the ODE in Eq. 15 is 
solved analytically by 

PR{t) = e''\p{0)+G-'K)-G-^K, (17) 
p{0) = Gs'K, 

where p(0) is the initial density matrix and e^* is the 
time-evolution superoperator for the density matrix. 

Calculating the exponential e^^ for a large number of 
time steps is computationally intensive, but we simplify 
by selecting a time-step resolution tgtep and using an it- 
erative process, 

Pi?(n*tstep) = (e^*^--)"(p(0) + G-iK)-G-^i^ . (18) 

One exponential is calculated for each selection of pa- 
rameters in G (including excitation frequency cj), and the 
problem is reduced to multiple matrix multiplications. In 
addition to calculating the matrix exponential, an inverse 

(These 
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matrix must also be calculated to solve Eq. 
two calculations prevent a general analytic solution and 
consume the most computational time.) The inverse of 
the steady-state superoperator Gs needs to be computed 
once for each selection of parameters, excluding the ex- 
citation frequency uj. 

Using these techniques we decrease the computation 
time of p{t) by three orders of magnitude compared to 
the conventional ODE solvers that were used in previous 
studieJ^^^"^. This makes the simulation of complex dis- 
tributions possible. For example, the data of Figs. [2] and |4] 
are superpositions of 2880 separate simulations generated 
at a resolution that would be impractical using conven- 
tional ODE solvers on a standard personal computer. We 
first verified the Liouville-space technique by successfully 
generating the uncoupled and exchange-cou pled s imula- 
tions previously generated using ODE solverPE^. Then 
the simulations obtained for dipolar-coupled pairs (see 
Figs, [l] and [3| are corroborated by ODE-based simula- 
tions (e.g., MATLAB® solver 0DE113). From our sim- 
ulations it is possible to describe the nature of the cou- 
pling within the pair that leads to experimentally ob- 
served spin-dependent transport and recombination pro- 
cesses. 



IV. RESULTS AND DISCUSSION 

The simulations are used to generate a representative 
database of different coupling strengths and Larmor sep- 
arations. Specifically, we discuss dipolar coupling within 
the intermediate-spin-pair model to account for \^jBi 
Rabi frequencies of experimental pODMR/pEDMR data 



in disordered semiconductors^^. However, we find 
that dipolar coupling alone does not account for certain 
data — exchange coupling must also be included. 

For each of the simulations a set of global parameters 
is used. Evolving p (r) during the application of a 2 /is 
excitation pulse, we calculate the observable Q{r) with 
a 4001-step resolution for a range of pulse frequencies cj, 
where the range of uj is covered with an 801-step resolu- 
tion. We choose ujo/27t = 10 GHz (within the microwave 
X-band) with the Larmor separation Acu centered on this 
value. For all simulations we also choose a Bi strength 
such that /27r = 10 MHz. The rate coefficients for 
singlet recombination, triplet recombination, and disso- 
ciation were assigned values of r^^ = 1 /is, r^^ = 100 /is, 
and = 1 ms, respectively. These rates represent the 
characteristic times for spin-dependent recombination or 
spin-independent dissociation into free charge carriers. 
To ensure all singlet information is recorded after the 
excitation, the observable (Eq. |4| is integrated up to a 
time to = 4r^^, with rs defined in Eq. pi This is done 
to offset the effects of the inherent signal reduction that 
arises as the exchange- or dipolar-coupling strength is in- 
creased. The generation rate k is chosen such that the 
initial (steady-state) pure-triplet populations of the den- 
sity matrix are approximately 0.05 (pfi(O) = P44(0) ^0.05 
in the 4x4 representation). 

All values are taken to be representative of measure- 
ment conditions that can be realized in the laboratory, 
following Ref. 18, Important physical information gar- 
nered from the simulations are the relative positions and 
amplitudes of the Rabi frequency components ft and their 
dependence on the different coupling strengths. 



A. Dipolar Coupling Only 

Here we vary the dipolar-coupling strength D with re- 
spect to the Larmor separation Auj^ with a negligibly 
small exchange interaction J. Displayed in Fig. [l] are 
simulations with Larmor separations of A(jj/27t = 1 MHz, 
20 MHz, and 40 MHz; mapped against dipolar coupling 
strengths i:>/27r = 1 MHz, 10 MHz, and 50 MHz. These 
are chosen as representative values of Aco'/27r and D/27r 
(smaller, approximately equal, and larger) relative to the 
excitation-field strength 75i/27r = 10 MHz. General fea- 
tures of these data include the resonance curves at 7^1, 
^/2JBl, and 27^1, which appear variously as a func- 
tion of Auj/27r and D/27T. The prominent vertical lines 
in Fig. [ijd) and (g) result from extremely long integra- 
tion times compounded with the continuous rotation into 
"leaky" singlet states. 

Weak dipolar coupling {D/2ti = 1 MHz) is shown in 
the top row of Fig. TJ(a)-(c)]; we reproduce qualitative 
features of the weakly coupled pair discussed in Ref. [171 
(The small differences from an uncoupled pair are the re- 
sult of a slight splitting of the resonances caused by weak 
dipolar coupling.) All plots in the top row of Fig. [l] have 
on-resonance single-transition Rabi frequencies of 7^1. 
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FIG. 1. (Color online) Plots of the Fast Fourier Transform FFT{Q(r)} of the observable Q{r) as a function of the excitation 
frequency cj, in the regime of dipolar coupling only. The signal intensity for each plot is normalized to plot (c) and given by the 
number next to the color scale, which indicates the highest magnitude signal intensity in the scale for that plot. Simulations 
are done with Larmor separations of Auj/27t = 1 MHz [plots (a), (d), and (g); left column], Auj/27t = 20 MHz [plots(b),(c), and 
(h); center column], and ALj/2n = 40 MHz [plots (c), (f), and (i); right column]; mapped against dipolar-coupling strengths of 
D/27T = 1 MHz [plots (a)-(c), first row], D/27T = 10 MHz [plots (d)-(f), second row], and D/h = 50 MHz [plots (g)-(i), third 
row]. The excitation strength is jBi/27i — 10 MHz. 



The multiple-transition Rabi frequencies arise from sim- 
ple addition and subtraction of the single-transition Rabi 
frequencies (see Ref. 17). Thus a weak-dipolar regime 
leads to no measurable V2jBi Rabi frequency compo- 
nents. Fig. [TJa) has an on-resonance Rabi frequency of 
27^1 due to a spin-beating effect from the coherent nuta- 
tion of both spin a^^ ' ^^ l Also, the intensity in Fig. TJa) is 
approximately half that of Fig. [TJc); this results directly 
from the relative triplet / singlet content of the eigenbasis 
in Eq. |3] The middle row of Fig. [l] has an intermediate- 
dipolar strength {D/2t: = 10 MHz = 7^1) and no strong 
\f2'^B\ Rabi frequency components. Indeed this ^/2jBi 
component is barely visible in Fig. [ijd), much weaker 
than the bright vertical lines. 



Strong dipolar coupling {D/27T = 50 MHz) is shown 
in the last row of Fig. [l] When the Larmor separa- 
tion is less than the dipolar strength, Auj < jBi < D 
[Fig. [ijg)], there is a weak non-visible (due to bin size) 
transition with a ^/2^Bl Rabi frequency. Both Fig. [l|h) 
and Fig. I]^i) show a strong ^/2JBl Rabi frequency com- 
ponent. Fig.jljh) has a dipolar-coupling strength greater 
than the Larmor separation, and both are greater than 
or comparable to the excitation strength, D > {Au ^ 
jBi). Fig. [TJi) has a Larmor frequency separation and 
dipolar-coupling strength approximately equal, but both 
greater than the excitation strength, {Auj ^ D) > jBi. 
Thus a Rabi frequency of V^7^i only occurs in the 
regime where dipolar-coupling strength is greater than 
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FIG. 2. (Color online) (a) The distribution of dipolar- 
coupling strengths for the simulation shown in (b) . The distri- 
bution is a Pake doublet with Larmor separation Auj/27t = 40 
MHz and dipolar-coupling strength of D/27r — 80 MHz con- 
voluted with a Lorentzian with a half-width of 10 MHz. (b) 
Plot of the Fast Fourier Transform FFT{(5(r)} of the observ- 
able Q(r) as a function of the excitation frequency uj. The 
signal intensity is normalized and given by the number next to 
the color scale, which indicates the highest magnitude signal 
intensity. The simulation uses the distribution in Fig. [2ja), 
with no exchange coupling J = 0. The excitation strength is 
7Bi/27r = 10 MHz. 



both the Larmor separation and the excitation strength, 
D > Auj^D > 7^1. The limits of this regime are dis- 
cussed in Sect. IIIIB2[ 

Each column of Fig. [l] reflects the observable intensity 
getting weaker with increasing dipolar-coupling strength; 
this is because we are approaching a triplet-singlet energy 
eigenbasis. Another trend occurring down each column 
is the separation of the on-resonance positions increasing 
with dipolar-coupling strength, also demonstrated with 
Eq.[lT] 

We turn now to modeling a more realistic 
pEDMR/pODMR signal. We operate in a regime 
of strong ^/2jBi Rabi components and seek to model 
an actual spin-pair distribution in a disordered material. 
Indeed, many materials with pronounced spin-selection 
rules are disordered semiconductors, including those for 
which the significance of the dipolar interaction has been 
discussed^^ In a disordered environment, the orien- 
tation of a spin pair with respect to an applied magnetic 
field can be entirely random. The strengths of the 



dipolar fields are highly orientation dependent because 
of the inherently anisotropic spin-dipolar interaction, 
even if a fixed spin-pair distance is considered (rather 
than a distribution of distances). The well-known Pake 
distribution accounts for this random orient at iorP^. 
Fig. |2] is a simulation using a Pake distribution with a 
dipolar coupling strength of I^/27r = 80 MHz, a Larmor 
separation of Acj/27r = 40 MHz, and J = 0. Fig. [2|a) is 
created using a 2880-point Pake distribution convolved 
with a 10 MHz FWHM Lorentzian function to account 
for power broadening due to the excitation pulse. We 
then generate simulations for the 2880 dipolar-coupling 
strengths and, using relative weights from Fig. [2|a), 
average those simulations to make Fig. [2|b). 

A comparison of Fig. [2|b) with experimental pEDMR 
and pODMR data^^ strongly supports the notion that 
the strong transitions with a Rabi frequency of a/275i 
arise from a strong dipolar interaction. However, other 
characteristics of Fig. |2|b) do not match experimen- 
tal data. The strong low-Rabi-frequency components 
OAjBi) of Fig. ^h) are not seen in pODMR of hy- 
drogenated amorphous silicon (a-Si:H), as seen in Fig. 8 
of Ref. 29 or Fig. la of Ref. |30l Low-Rabi-frequency 
components (~ 0.1-0.2 7^1) are seen in pEDMR of hy- 
drogenated amorphous silicon nitride (a-SiN^^iH), given 
in Fig. 2(e) of Ref. [311 However, the same data also 
shows a \/2'^B\ Rabi frequency relatively flat with re- 
spect to excitation frequency compared to the curved 
shape in Fig. [2|b). From these discrepancies we con- 
clude that dipolar coupling alone cannot account for the 
pODMR/pEDMR data reported in the literature. 



B. Dipolar and Exchange Coupling 

We now introduce exchange coupling between the spin 
pairs in addition to the dipolar coupling. The parameter 
space for the simulation is quickly growing; we give a 
small representation in Fig. [3j 

Fig. [3] has plots of the Fast Fourier Transform 
FFT{(5(r)} of the observable Q(t^ as a function of the 
excitation frequency uj. The signal intensity for each plot 
is normalized to plot (a) and given by the number next 
to the color scale, which indicates the highest magni- 
tude signal intensity in the scale for that plot. Simula- 
tions are done with a Larmor separation of Aa;/27r = 40 
MHz. Dipolar-coupling strengths are Dl2i\ = [plots 
Fig. [3];a),(eUi),(m),(q), first column], D/27r = 10 MHz 
[plots Fig. [3]^b),(f),(j),(n),(r), second column], D/2tt = 
40 MHz [plots Fig. [3];c),(g),(k),(o),(s), third column], 
D/2t: = 80 MHz [plots Fig. |];d),(h),(l),(p),(t), fourth 
column] ; mapped against exchange-coupling strengths of 
J/27r = [plots Fig. [3]; a)- (d), first row], J/27r = 10 MHz 
[plots Fig. [3];e)-(h), second row], J/27r = 50 MHz [plots 
Fig.[3];i)-(1), third row], J/2ir =80 MHz [plots Fig.Qm)- 
(p), fourth row], J/27r = 300 MHz [plots Fig. [3|;qj-(t), 
fifth row]. With the expection of Fig. [3]^ a), dipolar- and 
exchange-coupling strengths are chosen greater than or 
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FIG. 3. (Color online) See Sect. IV B for description. 




equal to the excitation strength 75i/27r = 10 MHz. 

Fig. [3]^ a) is an uncoupled spin pair that satisfies the 
weak-coupling limit described in Sect. |IIIB 1| two res- 
onances are located at the Larmor frequencies of the 
electron and hole and have Rabi frequencies of jBi. 
The uncoupled spin pair yields the maximum relative 
intensity (100) in Fig. [sj Fig. [sjf) and Fig. [sjp) are 
in the effectively-weak-coupling limit also described in 
Sect. |IIIB 1[ where the dipolar- and exchange-coupling 
strengths are equal {J = D). Fig.jsjk) has approximately 
equal dipolar- and exchange-coupling strengths with on- 
resonance Rabi frequencies slighty offset from jBi. 

The first row of Fig. [3] is similar to the third column of 
Fig. [1] there is no exchange interaction present and the 
relative intensity of FFT{(5(r)} decreases with increasing 
dipolar coupling strength. The distribution in Fig. [2|b) 
can be thought of as generated from intermediate val- 
ues between and including Fig. [3|(a)-(d)]. The sequence 
across the fir st row o f Fig. [3] best illustrates the discus- 



sion in Sect. IIIB2 



the two |T) ^ 1 2) transitions are 
split from the center frequency ujq and trend upwards to 
the strong dipolar-coupling limit with a Rabi fre- 

quency. The two |T)^|3) transitions are also split from 
the center frequency but are approaching their strong- 
coupling limit (zero Rabi frequency). 

As discussed in the strong-exchange-coupling limit of 
Sect. |IIIB 3[ the single-transition analysis fails to account 
for the observed Rabi frequencies; this is explicitly seen 
down the first column of Fig.jsj for the two |T) ^ |3) tran- 
sitions. Indeed the observed frequencies are not 
and as would be obtained from Eq. [lOj multiple tran- 
sitions must be considered to obtain the correct values. 
The first column in Fig. |3] is an extension of the simu- 
lations shown in the third column of Fig. 2 in Ref. [181 
where the exchange interaction is considered without the 
dipolar interaction. With increasing exchange-coupling 
strength the |T)^|2) transitions are split further about 
cjo, while the |T)^|3) on-resonance frequency positions 
remain unaffected. This is seen down each column of 
Fig. [3] and from Sects. [TTTb1| and |III B 4} The role dipo- 
lar coupling plays is shown by a common trend through- 
out all rows of Fig. [3) An increase of dipolar-coupling 
strength creates a greater energy splitting, causing a par- 
ticular transition to tend further from ujq and making the 
single-transition analysis of Eq. 10 valid. Therefore, a 
^/2jBi Rabi frequency is present if the dipolar coupling 
is strong enough. From these general trends, we deter- 
mine that only the combination of strong dipolar and 
even stronger exchange (see Fig. [3[(r)-(t)]) yields strong 
\^'yBi Rabi frequency components without any strong 
low-frequency (0-7^1) components. 

Using this analysis we generate Fig. [4]^b), which shows 
a distribution simulation similar to that of Fig. [2|b) but 
with a strong exchange coupling. This distribution sam- 
ples from the regime where there is a large difference 
between exchange- and dipolar-coupling strengths with 
J>D>Auj> jBl, J - D > Acj. Fig. gb) has a fiat 
Rabi frequency of yf^'^Bx and exhibits no strong low- 
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FIG. 4. (Color online) (a) The distribution of dipolar- 
coupling strengths for the simulation shown in (b) . The distri- 
bution is a Pake doublet with Larmor separation /\uj/2ti — 40 
MHz and dipolar-coupling strength of D/27r = 80 MHz convo- 
luted with a Lorentzian with a half-width of 10 MHz. (b) Plot 
of the Fast Fourier Transform FFT{(3(r)} of the observable 
Q(r) as a function of the excitation frequency uj. The signal 
intensity is normalized to Fig. [2|b) and given by the number 
next to the color scale, which indicates the highest magni- 
tude signal intensity. The simulation uses the distribution in 
Fig. [4|a), with an exchange-coupling strength of J/2ti — 300 
MHz. The excitation strength is 7Bi/27r = 10 MHz. 



Rabi- frequency components. It also exhibits some 27^1 
components. These same characteristics are found in the 
experimental data of Ref. [29ll3T 

The pODMR data of a-Si:H in Fig. 8 of Ref. |29lis al- 
most identical to Fig. [ij^b), showing strong ^/2^Bi Rabi- 
frequency components, weak components around 27^1, 
and no low-frequency components. Thus our spin-pair 
model predicts that both dipolar and exchange coupling 
are responsible for the pODMR data of Ref. ^29, More- 
over, the simulations show that the relative coupling 
strengths present in this data are in a regime with a large 
difference in exchange- and dipolar-coupling strengths, 
with J > D. This analysis supports the discussion pre- 
sented in Ref. ^ that suggested dipolar coupling was 
the cause for the observed data; it further predicts that 
strong exchange coupling was also present. 

The pODMR of a-Si:H geminate pairs in Fig. 3a of 
Ref. [3Q| is also very similar to the simulation in Fig.|4]^b), 
with the caveat that there appears to be the presence of 
uncoupled spins that produce strong jBi Rabi frequen- 
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cies. In that data set, the strong transitions with a Rabi 
frequency of \/2'^B\ are fiat with respect to excitation 
frequency but become abruptly weaker; this is character- 
istic of the Pake distribution in Fig. |4]^b) which also has 
a strong \/2^Bi component becoming abruptly weaker 
at an excitation frequency of {uj — ujq)/^Bi = 10. This 
experimental data also has no lower- Rabi-frequency com- 
ponents (0-7^1), which we have shown to be a defining 
characteristic of the regime in which there is a large dif- 
ference in exchange- and dipolar-coupling strengths with 
J>D. Therefore, we determine that the geminate pairs 
show the characteristics of uncoupled pairs mixed with 
strongly dipolar-coupled pairs discussed in Ref. 30, but 
we also predict the presence of a strong exchange cou- 
pling. 

Finally, the pEDMR data in Fig. 3e of Ref. ISTI show 
broad \/2'^Bi Rabi frequency components and weak 7^1 
Rabi frequency components. Again, this could be charac- 
teristic of a resonance involving mostly uncoupled pairs 
and some strong exchange- and dipolar-coupled pairs. 
However, the presence of both strong dipolar and ex- 
change coupling cannot explain the strong low-frequency 
0.27^1) components present in the Ref. 31 data. 
Perhaps the curvature leading to the \/2^Bi limit seen 
in Fig. [2] cannot be seen in the Ref. |3T] data because 
of a low number of dipolar-coupled pairs. However, if 
the strong low-frequency components are due to strong 
dipolar coupling alone, we would expect (from Fig. [2| 
that the y/2^Bi component would be as strong as the 
low-frequency component, and the data does not have 
this feature. Therefore we conclude that dipolar-coupled 
pairs can explain the y/2'yBi Rabi frequencies in the a- 
SiNa^iH data presented in Ref. 31, but whether exchange 
is present cannot be confirmed or rejected due to the 
weak signal strength of the strongly coupled pairs rela- 
tive to the uncoupled pairs in that data. 



V. SUMMARY AND CONCLUSION 

Numerical and analytical methods were used to in- 
vestigate the role of the dipolar interaction for electri- 
cally and optically detected Rabi oscillation frequencies 
of intermediate-spin-pair systems. A general descrip- 
tion of the physics of pEDMR and pODMR transient- 
nutation experiments was given that includes dipolar 
and exchange interactions, the Larmor separations within 
the intermediate pairs, and the excitation- field strength. 
We have presented a intermediate-spin-pair model that 
corroborates previous numerical studies that included 
weakly coupled pairs only^^ and exchange-coupled 
pairs onl}^. The model also supports experimental 
studies that attributed the observation of \/2^Bi Rabi- 
frequency components with pODMR/pEDMR of disor- 
dered semiconductor J^^M^ to the presence of strong dipo- 
lar coupling within the spin-pair model. We have shown 
that pODMR data of a-Si:H presented in Ref. 29 and [30] 
can be explained within a intermediate-pair model in the 
regime of strong dipolar coupling and stronger exchange 
coupling, (J - > Ao;^, J>D>Auj>'yBi. 
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